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I study the diffuse flux of electron antineutrinos from stellar collapses with direct black hole 
formation (failed supernovae). This flux is more energetic than that from successful supernovae, 
and therefore it might contribute substantially to the total diffuse flux above realistic detection 
thresholds. The total flux might be considerably higher than previously thought, and approach 
the sensitivity of SuperKamiokande. For more conservative values of the parameters, the flux from 
failed supernovae dominates for antineutrino energies above 30-45 MeV, with potential to give an 
observable spectral distortion at Megaton detectors. 



There is confidence, in the neutrino astrophysics com- 
munity, that the diffuse flux (DF) of neutrinos from 
core collapse supernovae will be detected in the near fu- 
ture. The current upper limit from the 50 kt water tank 
of SuperKamiokande on diffuse electron antineutrinos, 
<t> Ve {E > 19.3 MeV) < 1.4-2 cm^s" 1 at 90% confidence 
level [U [2] , already approaches theoretical predictions [3] . 
Within a decade or so, tens to hundreds of events from 
the diffuse flux will be available from detectors of 0.1-1 
Mt mass [4j [5] , allowing steady progress in the investiga- 
tion of the physics of core collapse, of the cosmological 
rate of supernovae and of the properties of the neutrino. 

So far, predictions for the DF have considered only the 
most common scenario: the collapse into a neutron star, 
with ~ 3-10 53 ergs emitted in neutrinos of average energy 
E ~ 9 - 18 MeV (see Eq. p)). Recently, detailed stud- 



ies have appeared [6j [7] 



, [TOl [TT] on the rarer case 



of direct collapse into a black hole without explosion, 
i.e., a failed supernova. It was shown that the neutrino 
emission is somewhat more luminous and decidedly more 
energetic than for neutron star- forming collapse, due to 
the rapid contraction of the newly formed protoneutron 
star preceding the black hole formation. Average energies 
are Eq ~ 20 — 24 MeV for all neutrino flavors. This sug- 
gests that the hotter contribution of black hole-forming 
collapses to the DF might exceed that of neutron star- 
forming ones in part of the energy spectrum. 

In this letter I study the diffuse flux from failed super- 
novae, with focus on its the electron antineutrino com- 
ponent, which is relevant for water Cherenkov detectors. 
The results confirm the intuition that the flux from failed 
supernovae might be significant and bring the DF even 
closer to being finally detected. 

Core collapse occurs for stars with mass M ^ 8M 
(Mq is the mass of the Sun), at an average rate of 
R cc (0) ~ 10~ 4 Mpc" 3 yr- 1 today [12 and of 

R cc {z)~R cc (0)(l + zf (1) 

(best fit (3 = 3.28 [12 J, up to redshift z ~ 1. The rate 
flattens at larger z. For M = 8 — 25M [13] the collapse 
leads to an explosion - the observed explosion of ~2OM 
star Sanduleak into SN1987A [14 supports this - and to 
the formation of a neutron star. 



Since stars are distributed in mass as (j)(M) ex M -2 35 
[15], one gets that neutron star- forming collapses are a 
fraction f^s — 0-78 of the total. Their neutrino output 
is roughly equipartitioned between the six species: v ei v e 
? ^/i > ^/x ? ? (^5^'^'^ = ^ from here on). At the 
production site, the flux in each species w, differential in 
energy, can be described as [16] : 
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where T(x) stands for the Gamma function. Here a w 
controls the spectral shape, L w is the time integrated 
luminosity and Eq w is the average energy. For illustration 
(fig. [I] dashed lines), here I use the typical values [T6j : 
E e = 15 MeV, E 0x = 18 MeV, L E = L x = 5 • 10 52 ergs, 
«e = 3.5 and a x = 2.5. After neutrino oscillations in the 
star (those in the Earth are negligible [3 ), the v e flux is 
determined by the survival probability p: 

F- e = pF9 + (1 - p)F x J , p = 0- cos 2 6>i2 ~ - 0.68 ,(3) 

with #i2 — 34° [17]. The interval for p reflects the possi- 
ble different oscillation scenarios, depending on the neu- 
trino mass hierarchy, on the still unknown angle #13 [18] 
and on the effects of neutrino-neutrino scattering (see 
e.g., 033 [20]). p can change during neutrino emission 
[2T] , but always remains in the interval Q , which there- 
fore describes the time averaged probability. To illustrate 
how results can vary, I consider the extreme values of p. 
They can be taken as constant in energy, which is ade- 
quate [18] for the experimentally relevant energy range 
(E > 19.3 MeV for SuperKamiokande ffl). 

While neutron star-forming collapses (NSFCs) have 
been widely studied, the evolution of higher mass stars is 
more uncertain. For M ~ 25-40 M (13% of the total) a 
weaker explosion should occur, with a black hole formed 
by fallback [13]. Stars with M £ 40 M (a 9% frac- 
tion), would instead collapse into a black hole directly. 
Simulations of such direct black hole-forming collapses 
(DBHFCs) H [9] [10] show an emitted neutrino flux 
that is more energetic and more luminous than the NSFC 
case, with especially high luminosity in v e and v e due to 
capture of electrons and positrons on nucleons. 
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FIG. 1: Neutrino fluxes at production inside the star for direct black hole- forming collapse (solid, from [TT]) , and neutron 
star-forming collapse (dashed, Eq. |2])). In both cases, the curves from upper to lower at 5 MeV correspond to zy e , z> e , v x . For 
direct black hole-forming collapse the neutrino spectra are shown for the Shen et al. (left panel) and Lattimer-Swesty (right) 
EoS. For each, the neutrino luminosities and average energies are given (inserts). See text for details. 
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A "stiffer" equation of state (EoS) of nuclear matter [7] 
and/or a smaller accretion rate of matter on the pro- 
toneutron star [9] correspond to more luminous and hot- 
ter neutrinos. Here I use the fluxes from DBHFCs as in 
fig. 5 of ref. [11 . They are are shown in fig. [I] (solid 
lines and inserts). These fluxes were obtained for the 
4OM progenitor in [22] with the stiffer Shen et al. (S) 
EoS [23] (incompressibility K = 281 MeV) and the softer 
Lattimer-Swesty (LS) one [24 (with K = 180 MeV [TT]). 
For the different progenitors considered in [TT] results ap- 
pear unchanged for the S EoS, while for the LS one the 
luminosity and average energy may be lower by a factor of 
two and by 10-20% respectively. For the energy spectra, 
I use the same linear interpolation of numerically calcu- 



lated points as in [11 , which underestimates the DF in 
the SuperKamiokande window by about 10-20%, so my 
results are conservative. Ref. [11 shows that Eq. |3| ap- 
plies to the DBHFC case, with the same extreme values 
of p being realized for the same oscillation parameters as 
in the NSFC case. 

I calculated the neutrino fluxes from NSFCs and DB- 
HFCs, and the total DF for a schematic two-population 
scenario, with a fraction f^s (Ibh — 1 — Ins) of iden- 
tical neutrino emitters of the NSFC (DBHFC) type [31]. 
Generalizing the single population formula (e.g., [25 ) one 
gets the total diffuse v e flux at Earth, differential in en- 
ergy and area: 



$(E) = — / R cc {z) [f NS F s NS (E(l + z)) + (1 - f NS )F s BH (E(l + z))} 
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where Q m = 0.3 and = 0.7 are the fractions of the 
cosmic energy density in matter and dark energy; c is 
the speed of light and Hq is the Hubble constant. I took 
the parameters R cc (0) = 10~ 4 Mpc- 3 yr-\ (3 = 3.28 
and z max = 4.5 [12 (results depend weakly on z maxi 
at the level of ~ 7% or less for z max > 3 [25 J. To 
parametrize the uncertainty in f^s I take the interval 
Ins = 0.78 — 0.91, corresponding to a mass 25 — 40 M 
as upper limit for neutron star- forming collapse. 

Results are shown in fig. [2] The diffuse flux from NS- 
FCs, ^ats, is maximum at E ~ 5 — 6 MeV, with an ex- 



ponential decay at higher energy [26]. The contribution 
from DBHFCs, has hotter spectrum, and thus is 

increasingly important at higher energy. Oppositely to 
*&NSi ^bh is larger for minimal permutation (p = 0.68) 
[TT] , because of the high original v e flux. The dependence 
of the original fluxes on the EoS is evident in $>bh- 

Fig. [2] shows that &bh might dominate already at E ~ 
22 MeV, implying a strong effect at SuperKamiokande. 
For the most favorable parameters the total flux above 
19.3 MeV more than doubles compared to 100% NSFCs, 
(fig. [3|, reaching a value ($ ~ 0.89 cm _2 s _1 ) tantaliz- 
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ingly close to the current upper limit. The enhancement 
of the event rate is even larger, thanks to the ~ E 2 depen- 
dence of the detection cross section. Thus, the DF might 
be closer to detection than previously thought, within 
the reach of improved searches at SuperKamiokande. 

It is more likely, however, that &bh exceeds $ns only 
above 30-40 MeV (fig. |2j. Its effect would be below the 
sensitivity of SuperKamiokande - which would therefore 
place limits on neutrinos from DBHFCs - but might be 



visible with the 1 Mt planned Cherenkov detectors, where 
(9(10) events are expected in this energy interval for a 
few years running time. Besides the excess in event rate, 
which suffers normalization uncertainties, the DBHFC 
diffuse flux could be visible for the spectral distortion 
that it produces. The lower threshold (~ 11 MeV) of a 
liquid scintillator [27 or Gadolinium-loaded water detec- 
tor [28] could allow to see a break in the energy spectrum 
at ~20 MeV, that might escape a pure water detector. 
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FIG. 2: The diffuse flux of v e at Earth from neutron star- forming collapses (dashed curves) and direct black hole- forming ones 
(solid), for different EoS and survival probability p. Direct black hole -forming collapses are assumed to be 22% (thick curves) 
or 9% (thin curves) of the total (see Eq. Q). 



If all parameters conspire to maximally suppress it, 
$ bh might be invisible even for a Mt detector at least in 
the first few years running time. A negative result would 
then constrain the parameter space strongly. 

To illustrate how results change for a more energetic 
NSFC flux I have repeated the calculations with Eq x = 22 
MeV; in what follows I discuss how results compare to 
those in fig. [2] For p — 0.68 differences are only minimal 
relative to fig. [2] This is because the v e flux at Earth is 
dominated by the original v e flux, which for the DBHFCs 
is markedly more energetic than in the NSFC case. For 
the S EoS <$>bh exceeds <$> NS above 22 MeV (38 MeV) 
for f NS = 0.78 (f NS = 0.91). Instead, for the LS EoS 
the two components are comparable at energy of 44 MeV 
or higher. This slight worsening compared to fig. [2] may 



make the difference between a positive or negative signal 
of &bh at an experiment. For p = the distinctive 
original v e flux from DBHFC does not contribute to the 
v e flux at Earth, and so the main signature of direct black 
hole formation is lost. One may see a slight excess flux 
at E ~ 45 - 50 MeV only for the S EoS and f NS = 0.78. 

If the flux excess due to &bh appears only above 30 
MeV, it might partially be masked by the invisible muon 
and atmospheric neutrino backgrounds, which are strong 
at that energy pQ. I modeled these following [29 and 
taking a 100% flux excess in the 30-35 MeV bin, relative 
to a theoretical model or fit to the data with f^s = 1- 
For pure water, where invisible muons dominate, the ex- 
cess would not be statistically significant in the single 
bin, but might be distinguishable with a fit of the spec- 
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trum of events. With the reduced background allowed 
by Gadolinium, 3a significance in the single bin would 
be achieved with about 12 Mt-yr exposure, with a lower 
exposure needed if a spectral fit is done. 

To summarize, the diffuse flux of neutrinos from failed 
supernovae may be significant, at a level detectable at 
SuperKamiokande or at Mt scale detectors. While con- 



clusions are limited by uncertainties, it is hoped that this 
letter will serve as a motivation for the development of 
more realistic predictions of the diffuse flux, which would 
be of great service to the experimental community. 

I thank H.T. Janka, A. Mezzacappa and O. L. G. Peres 
for useful discussions, and ASU and RBRC for support. 
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FIG. 3: Left: the largest possible v e diffuse flux from direct black hole-forming collapses (solid, thick line), compared to that 
from neutron star-forming ones (dashed), from fig. [5] The total flux is shown too (thin). Integrated fluxes above two thresholds 
of interest are given, in cm _2 s _1 . Right: same figure for inverse beta decay events in water with 2 Mt-yr exposure. E e is the 
positron energy. I used R cc (0) = 10 -4 Mpc _3 yr _1 for the core collapse rate today and ref. [30] for the detection cross section. 
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